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An effect of the differential rotation on the nonlinear electromotive force in MHD turbulence is 
found. It includes a nonhelical a effect which is caused by a differential rotation, and it is independent 
of a hydrodynamic helicity. There is no quenching of this effect contrary to the quenching of 
the usual a effect caused by a hydrodynamic helicity. The nonhelical a effect vanishes when the 
rotation is constant on the cylinders which are parallel to the rotation axis. The mean differential 
rotation creates also the shear-current effect which changes its sign with the nonlinear growth 
of the mean magnetic field. However, there is no quenching of this effect. These phenomena 
determine the nonlinear evolution of the mean magnetic field. An effect of a uniform rotation on 
the nonlinear electromotive force is also studied. A nonlinear theory of the fix J effect is developed, 
and the quenching of the hydrodynamic part of the usual a effect which is caused by a uniform 
rotation and inhomogeneity of turbulence, is found. Other contributions of a uniform rotation to 
the nonlinear electromotive force are also determined. All these effects are studied using the spectral 
t approximation (the third-order closure procedure). An axisymmetric mean-field dynamo in the 
spherical and cylindrical geometries is considered. The nonlinear saturation mechanism based on 
the magnetic helicity evolution is discussed. It is shown that this universal mechanism is nearly 
independent of the form of the flux of magnetic helicity, and it requires only a nonzero flux of 
magnetic helicity. Astrophysical applications of these effects are discussed. 

PACS numbers: 47.65.+a; 47.27.-i 



I. INTRODUCTION 

Generation of magnetic fields by a turbulent flow of 
conducting fluid is a fundamental problem which has 
a large number of applications in solar physics, astro- 
physics, geophysics-planetary physics and in laboratory 
studies (see, e.g., |1 H 1 H BE II II @, and refer- 
ences therein). In recent time the problem of nonlinear 
mean-field magnetic d yna mo is a subject of active dis- 
cussions (seeje.g;, PS 111 EEQIllEmEE 
l2fl l2l l2l l2l |24|. and references therein). The con- 
ventional approach to the nonlinear dynamo is based on 
comparison of the three effects participating in dynamo 
action, namely the a effect (caused by helical motions 
of a turbulent fluid) , the large-scale differential (nonuni- 
form) rotation <5f2 and the turbulent magnetic diffusivity 
rj T . The mean magnetic field is generated due to a com- 
bined effect of the differential rotation and the a effect. 
These effects have been considered as independent phe- 
nomena. In particular, the electromotive force has been 
determined independently of the differential rotation. 

On the other hand, the differential rotation can be re- 
garded as large-scale motions with a mean velocity shear 
imposed on the small-scale turbulent fluid flow. An in- 
teraction of the mean differential rotation with the small- 
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scale turbulent motions can cause a generation of a mean 
magnetic field even in a nonhelical, homogeneous and 
incompressible turbulent fluid flow. This mechanism of 
mean-field dynamo is associated with a shear-current ef- 
fect which is determined by the WxJ term in the elec- 
tromotive force, where W is the mean vorticity caused by 
the mean velocity shear and J is the mean electric current 
(see [25| ). A nonlinear theory of a shear-current effect 
in a nonrotating homogeneous and nonhelical turbulence 
with an imposed mean velocity shear in a plane geome- 
try was developed in pfij . It was shown that during the 
nonlinear growth of the mean magnetic field, the shear- 
current effect changes its sign, but there is no quenching 
of this effect contrary to the quenching of the usual a 
effect, the nonlinear turbulent magnetic diffusion, etc. 

In this study we investigated the effects of differential 
and uniform rotation on nonlinear electromotive force. 
The main conclusion of this study is that the nonlinear 
electromotive force cannot be determined independently 
of the mean differential rotation. We found a nonheli- 
cal a effect which is caused by a differential rotation and 
is independent of a hydrodynamic helicity. There is no 
quenching of this effect contrary to the quenching of the 
usual a effect caused by a hydrodynamic helicity. The 
mean differential rotation of fluid can decrease the to- 
tal a effect due to the nonhelical a effect. Two kinds of 
the a effect (helical and nonhelical) have opposite signs. 
Therefore, the total a effect should always change its sign 
during the nonlinear growth of the mean magnetic field 
because there is a quenching of the usual (helical) a ef- 
fect. This can saturate the growth of the mean magnetic 
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field. 

The mean differential rotation creates also the shear- 
current effect. We found that the mean differential ro- 
tation increases the growth rate of the large-scale dy- 
namo instability at a weak mean magnetic field due to 
the shear-current effect, and causes a saturation of the 
growth of the mean magnetic field at a stronger field. 
Note that the applications of the obtained results to the 
solar convective zone shows that the nonlinear shear- 
current effect becomes dominant at least at the base of 
the convective zone. We found that the nonlinear func- 
tion (7 (B) defining the shear-current effect is the same 
for a turbulence with a mean differential rotation in cylin- 
drical and spherical geometries for an axisymmetric mean 
field dynamo problem and for a nonrotating turbulence 
with an imposed linear mean velocity shear in a plane 
geometry. The latter case was investigated in 26] . 

We also studied an effect of a uniform rotation on the 
nonlinear electromotive force. In particular, we devel- 
oped a nonlinear theory of the fix J effect and we de- 
termined the nonlinear quenching of the hydrodynamic 
part of the a effect which is caused by both, a uniform 
rotation and inhomogeneity of turbulence. Other nonlin- 
ear coefficients defining the nonlinear electromotive force 
are also determined as a function of a uniform rotation. 
In this study we considered a uniform rotation with a 
small rotation rate in comparison with the correlation 
time of the fluid turbulent velocity field. We studied all 
the above effects using the spectral r approximation (the 
third-order closure procedure). 

This paper is organized as follows. In Section II we for- 
mulated the assumptions and the method of the deriva- 
tion of the nonlinear electromotive force in a turbulence 
with a uniform and nonuniform rotations. In Section 
III we considered axisymmetric mean-field dynamo equa- 
tions and determined the coefficients defining the electro- 
motive force for a rotating turbulence. In Section III we 
also discussed in details the effects of differential and uni- 
form rotation on nonlinear coefficient defining the elec- 
tromotive force. In Section IV we analyzed the nonlinear 
saturation of the mean magnetic field and discussed the 
astrophysical applications of the obtained results. In Ap- 
pendix A we derived the nonlinear electromotive force in 
a turbulence with uniform and nonuniform rotations. 



II. THE METHOD OF DERIVATIONS 

In a framework of the mean-field approach the evo- 
lution of the mean magnetic field B is determined by 
equation 



of the electromotive force £ = (u x b) in an anisotropic 
turbulence is given by 



— = V x (U x B 

dt x 



S - r;V x B) 



(1) 



(see, e.g., |l|, |2j, y, U, |3, |6(), where U is a mean velocity 
(the differential rotation) , r\ is the magnetic diffusion due 
to the electrical conductivity of fluid. The general form 



£i = a l] B 3 + (V ett xB) l -7 llJ (VxB) J 
-[5x(VxB)]j - K ljk (dB) jk 



(2) 



(see [27|, |2g), where (dB) i3 = (1/2) (V^- + V,-^), u 
and b are fluctuations of the velocity and magnetic field, 
respectively, angular brackets denote averaging over an 
ensemble of turbulent fluctuations, the tensors and 
rjij describe the a-effect and the turbulent magnetic dif- 
fusion, respectively, V is the effective diamagnetic (or 
paramagnetic) velocity, Kijk and S describe an evolution 
of the mean magnetic field in an anisotropic turbulence. 
Nonlinearities in the mean-field dynamo imply depen- 
dencies of the coefficients (a,-j, rjij, V eff , etc.) defining 
the electromotive force on the mean magnetic field. 

The method of the derivation of equation for the non- 
linear electromotive force in a rotating turbulence is sim- 
ilar to that used in [26] for a nonrotating turbulence with 
an imposed mean velocity shear. We consider the case 
of large hydrodynamic and magnetic Reynolds numbers. 
The momentum equation and the induction equation for 
the turbulent fields in a frame rotating with an angular 
velocity ft are given by 



du(t, x) 
dt 



db(f,x) 
dt 



_ Vftot + J_[( b . V )B + (B ■ V)b] 

Po MPo 
+2uxft - (U • V)u - (u • V)U 
+u N + F , (3) 

(B • V)u - (u • V)B - (U • V)b 

+ (b-V)U + b Ar , (4) 



where V • u = 0, po is the fluid density, p is the mag- 
netic permeability of the fluid, pqF is a random exter- 
nal stirring force, u N and h N are the nonlinear terms 
which include the molecular dissipative terms, ptot = 
p + p~ Y (B • b) are fluctuations of the total pressure, p 
are fluctuations of the fluid pressure. Hereafter we omit 
the magnetic permeability of the fluid, p, in equations, 
i.e., we include /i -1 / 2 in the definition of magnetic field. 
We study the effect of a mean rotation of the fluid on 
the nonlinear electromotive force. We split rotation into 
uniform and differential parts. By means of Eqs. ©-0 
written in a Fourier space we derive equations for the cor- 
relation functions of the velocity field /y(k) =!/(«,; Uj), 
of the magnetic field /iy(k) = L(bi\ bj) and for the cross 
helicity <faj-(k) = L(bi\Uj), where 



L(a; c)=J (o(t, k + K/2)c(f, -k + K/2)) 

x cxp (iK-R) dK , (5) 

and R and K correspond to the larg e scales, and r and k 
to the small ones (see, e.g., pjj l30j ]). The equations for 
these correlation functions are given by Eqs. l|AT )l -<|A3 )l 
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in Appendix A. These equations for the second moments 
contain high moments and a closure problem arises (see, 
e -g-i [31Ll32 . l33j ]). We apply the spectral r approximation 
or the third-order closure procedure (see, e.g., [T^. I20L I3TL 
0jESEEI13)I13)) which allows to express the deviations 
of the third moments from the background turbulence in 
k space in terms of the corresponding deviations of the 
second moments, e.g., 



Df N 



Df„ 



N(0) 



JV(O) _ 



-(/«-4 0) )/r(fc), 
-9ij/T(k) , 



(6) 
(7) 
(8) 



where the tensors Dflj , DhM and Dgfj are related to 
the third moments in equations for the second moments 
fij,hij and gij, respectively (see Eqs. (|A1|) - (|A3|) in Ap- 
pendix A). The correlation functions with the superscript 
(0) determine the background turbulence (with a zero 
mean magnetic field, B = 0), and h[j is the nonheli- 
cal part of the tensor of magnetic fluctuations of the 
background turbulence, r(fc) is the characteristic relax- 
ation time of the statistical moments. We applied the 
r-approximation only for the nonhclical part hij of the 

tensor of magnetic fluctuations. The helical part 
depends on the magnetic helicity, and it is determined 
by the dynamic equation which follows from the mag- 
netic helicityconservation arguments 0, |3j| (see also 
[HQ El El HI E3) In the present paper we con- 
sider an intermediate nonlinearity which implies that the 
mean magnetic field is not enough strong in order to af- 
fect the correlation time of turbulent velocity field. We 
also consider uniform rotation with a small rotation rate 
in comparison with the correlation time of the fluid tur- 
bulent velocity field (SI To C 1). The mean velocity shear 
due to the differential rotation is considered to be weak 
(SUto <C 1). For the integration in k-space of the second 
moments we use the following model of the background 
turbulence (with zero mean magnetic field, B = and 
without rotation): 



4 0) (k) = E(k)i (u 



, 2 \(°) 



5ij kij + (fcjAj 



(v) 



1 



p 



Coo 



= (b 2 )^E(k) 



(9) 



(10) 



where e^-fe is the Levi-Civita tensor, <5y is the Kro- 
necker tensor, kij = kikj/k 2 , E(k) = ~(d?(k)/dk)/8nk 2 , 
r(k) = 2r f(fc), f(k) = (k/ko) 1 ^, 1 < q < 3 is the 
exponent of the kinetic energy spectrum (e.g., q = 5/3 



for Kolmogorov spectrum), fco = 1/Zo, and lo is the 
maximum scale of turbulent motions, r = Iq/u , uq 
is the characteristic turbulent velocity in the scale lo, 



Af J = V,(u 2 )(°)/(u 2 )(°), A? } = V^b 2 )^/^ 2 )^, and 
(u • (Vxu))(°) is the hydrodynamic helicity of 
the background turbulence, / f^(k) dk = ((u 2 )(°) /3)6 tj 



and J h^ 1 (k) dk = ((b 2 )(°)/3)% Note that 5 J } (k) = 0. 
Here we neglected a very small magnetic helicity in the 
background turbulence. However, the magnetic helic- 
ity in a turbulence with a nonzero mean magnetic field 
is not small (see Section III-D). The derived equations 
allow us to determine the nonlinear electromotive force 
£i = £imn J 9nm(k) dk in a rotating turbulence (see for 
details, Appendix A). 



III. THE NONLINEAR ELECTROMOTIVE 
FORCE IN A ROTATING TURBULENCE FOR 
AN AXISYMMETRIC DYNAMO 

We consider the axisymmetric aJl-dynamo problem. 
In cylindrical coordinates (p, tp, z) the axisymmetric 
mean magnetic field, B = B(p,z)e ip + V x[A(p, z)e v ], 
is determined by the dimensionless equations 



dA 
~dt 



dB 
~dt 



a(B) +W, 01(B) V,(«i) 



B + Va (B)A s A 



-(V A (B) ■ V)(pA) - W, 00(B) (ClB) 
P 

f2* 5%(B) (Cj ■ V)B , 



D (flA) + pV 



1 



fo B V-V B CB)](pB) 



(11) 



(12) 



where SI, = 3nr Q /R a = [fi/(<5fi)„] , W* = 
(h/L) 2 (R u /R a ), & = n/f2and 



(ClB) = [V z (6Sl)V p - V p (SSl)V z ](pB) , 
Va(B) 



MB) 



V B (B)=V d (B) + M * ) + M * ) e P , 

P 

V d (B) = -M^( A W- £ aW)A 



and A s = A-l/p 2 , and = (VB 2 )/B 2 . The nonlin- 
ear coefficients a(B), t) a (B), T] B (B) defining the nonlin- 
ear a effect and the nonlinear turbulent magnetic diffu- 
sion of the poloidal and toroidal components of the mean 
magnetic field, are determined by Eqs. Ijl5(l and (|21|) in 
Section III-A. The nonlinear coefficients <tq(B) defining 
the shear-current effect and cr\(B) defining the nonhcli- 
cal a effect, are determined in Section III-B. The coeffi- 
cient Sq(B) defining the nonlinear fix J effect, is deter- 
mined in Section III-C. The quenching functions <f> n (B) 
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are determined by Eqs. (1221) in Section III-A. Note that 
in the equations for the nonlinear effective drift velocities 
Vyi(B) and Vb(B) of the poloidal and toroidal compo- 
nents of the mean magnetic field we neglected small con- 
tributions ~ O[(lo/L) 2 ] caused by the mean differential 
rotation. 

Equations and 1|12|) are written in the dimen- 

sionless form, where length is measured in units of L, 
time in units of L 2 jr\ T and the mean magnetic field B 
is measured in units of the equipartition energy B eq = 
^/PoUq, the magnetic potential A is measured in units 



of R a L B eq , the nonlinear a is measured in units of a* 
(the maximum value of the hydrodynamic part of the 
a effect), the basic scale of the turbulent motions I and 
turbulent velocity •*/ (u 2 ) at the scale / are measured in 
units of their maximum values Iq and uq, respectively, the 
dimensionless parameters A^) and are measured in 
the units of L~ x and A( s ) is measured in the units of 
L , the differential rotation 5fl is measured in units 
of (<5f2)*, the nonlinear turbulent magnetic diffusion co- 
efficients r] A B (B) are measured in the units of rj T and 
the nonlinear effective drift velocities ~Va,b(B) are mea- 
sured in the units of r\ T jL. We define R a — La*/rj T1 
Ru; = (<5f2)* L 2 /r/ T , the characteristic value of the tur- 
bulent magnetic diffusivity r\ T = IqUq/3, the dynamo 
number D = R w R a and Rm = loUo/rj is the magnetic 
Reynolds number. 

In spherical coordinates (r, 9, ip) the axisymmetric 
mean magnetic field, B = B(r,8)e v + Vx [A(r, 9)e v }, 
is determined by the dimensionless equations 



OA 
~dt 



dB 
~dt 



a(B) +W,*i (B) V, (Sfl) 
1 



B + Va (B)A s A 
(V A (B)-V)A-W*a (B)(nB) 



-fil*«5?(B) (u>- V)B, 
D(QA) + r sinflV 

-V B (B)]B , 



where A — r sin 9 A, B — r sin 9 B, 



r 2 sin 2 



(13) 



(14) 



(riB) 

Va(B) 
V b (B) 



[V r (50) V e - V e (<5fi) V r ]B , 
cos 9 V r — s'm9We , 

V d (B) - Mil A w _ M^l 



(e r + cot 9 eg), 



A s = A— l/(r sin(9) 2 and V e = (1/r) (<9/<90). Note that 
p = r sin#. 



A. The nonlinear a effect and the nonlinear 
turbulent magnetic diffusion coefficients of the mean 
magnetic field 

The nonlinear a effect is given by a(B) = a v + a m , 
where a" — x V( t )V {B) + aU 1S the hydrodynamic part 
of the a effect, and a m = x c (B)0 m (-B) is the magnetic 
part of the a effect, and the dimensionless parameter 
X v = —t [i v /3a* is related to the hydrodynamic helicity 
/j 11 = (u ■ (Vxu))' ' of the background turbulence, the 
dimensionless function x c (B) = (r /3poa*)(b • (Vxb)) 
is related to the current helicity (b • (Vxb)). Here \ v 
and x° are measured in units of a*, tq — Iq/uq is the 
correlation time of turbulent velocity field and a Q is the 
contribution to the hydrodynamic part of the a effect 
caused by a uniform rotation and inhomogeneity of tur- 
bulence. Thus, 



a(B) = xVCB) + a n + x c (B)<T(£) 



[see Eqs. 



a 



28J) and 

2 L £7* 

3 1,, 



(15) 

50j) in Appendix A], where 
[^(fl)AW+ e ^(B)A(»)] ,(16) 



the quenching functions (j) v (B) and cj> m (B) are given by 



12 

7/32 

_3_ 

w 



1 



arctan (3 
arctan (3 



(17) 
(18) 



(see 0), where/3 = and L(y)_= 1 
y" 2 ). Thus = 2 1(3 2 and m (£) = 

and 4> V {B) = 1 - (6/5)/3 2 and <j> m (B) = 



-2y 2 + 2y 4 ln(l + 
3//3 2 for j8 > 1; 
1 - (3/5)/? 2 for 



/3 <C 1. The quenching functions 4>i{B) and (/)§(B) are 
determined by Eqs. (|A51|) and (|A52I) in Appendix A. 
The function x c (B) entering the magnetic part of the 
a effect is determined by the dynamical equation i|29|l . 
Note that in Eq. JTJjJ we neglected small contributions 
~ 0(Sft/Cl) caused by the mean differential rotation and 
inhomogeneity of turbulence [these effects are given by 
Eqs. (| A68|> - <| A70fl in Appendix A]. For a nonhelical back- 
ground turbulence the first term, x V( t> v {B), in Eq. 115|) 
vanishes. 

The contribution to the nonlinear a effect caused by 
a uniform rotation for a weak mean magnetic field B <C 
B oq /4 is given by 



16 LCI 



15 L n 



aW--a(') 

3 



3e 



(19) 



and B ^> -B oq /4 it is given by 



1 LCt, 



3/3 2 L 7 



w-(A 



:AW) 



^(^))(l-l. 3e) , 



(20) 
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[see Eqs. (|A53|) and l|A56|) in Appendix A], where the 
parameter e = (b 2 )(°'/(u 2 )(°) is the ratio of the mag- 
netic and kinetic energies in the background turbulence. 
Asymptotic formula (|19fl for a n in the limit of a very 
small mean magnetic field coincides with that obtained 
in 28] for q = 5/3. 

The splitting of the nonlinear a effect into the hydrody- 
namic, a v , and magnetic, a m , parts was first suggested in 
[34| . The magnetic part a m includes two types of nonlin- 
earity: the algebraic quenching described by the function 
(f> m (B) (see an d the dynamic nonlinearity which 

is determined by Eq. (|29|l . The algebraic quenching of 
the a-cffcct is caused by the direct and indirect modifi- 
cation of the electromotive force by the mean magnetic 
field. The indirect modification of the electromotive force 
is caused by the effect of the mean magnetic field on the 
velocity fluctuations and on the magnetic fluctuations, 
while the direct modification is due to the effect of the 
mean magnetic field on the cross-helicity (see [l9l |20| V 

The nonlinear turbulent magnetic diffusion coefficients 
of the mean magnetic field are given by 

Va (B) = MB), V B (B) = MB) + MB) (21) 

(see |2(|), where the quenching functions (f)k(B) are given 
by 

MB) = 4 1) (4S) + 4 1) (4S), (22) 
MB) = ~(l + e)A^(4B), 

MB) = (2-3e)4 1) (4B)-(l-e)|-A 2 (16B 2 ), 

Z7T 

the functions At(y) and A^\y) are given by Eqs. I|A22|) - 
(|A24|) in Appendix A. The asymptotic formulas for the 
functions (f>k(B) for B <C i3 eq /4 are given by 4>i(B) — 
1 - (12/5) /3 2 , MB) = 1 - (4/5) (1 + e) (3 2 and MB) = 
-(8/5) (1 - 2e) 1 . For B > B eq /4 they are given by 
MB) = I//? 2 , MB) = MB) = 2(1 + e)/3/3, where 
(3 = \/8B. Note that in Eq. l|21|l we neglected small con- 
tributions ~ O[(lo/L) 2 ) caused by the mean differential 
rotation. 



B. The nonlinear coefficients ao(B) and <Ji(B) 
denning the shear-current effect and the nonhelical 
a effect 

The nonlinear coefficient <Jq(B) describes the shear- 
current effect (see [2^, |2^) and 0i(B) determines the 
nonhelical a effect. The parameters uq(B) and pi(B) are 
determined by the corresponding contributions from the 
J(B) term, the rjij (B) term and the Kyfc(B) term in the 
nonlinear electromotive force J5J) caused by the mean dif- 
ferential rotation. We found that the nonlinear function 
(7q(B) defining the shear-current effect is the same for a 
turbulence with a mean differential rotation in cylindri- 
cal and spherical geometries for an axisymmetric mean 
field dynamo problem and for a nonrotating turbulence 



with an imposed linear mean velocity shear in a plane 
geometry. The latter case was studied in 26] . 

To explain the physics of the shear-current effect, we 
compare the a effect in the aO. dynamo with the 5 term 
caused by the shear-current effect (see [25J |2(|). The 
a term in the nonlinear electromotive force which is re- 
sponsible for the generation of the mean magnetic field 
and caused by a uniform rotation and inhomogeneity of 
turbulence, reads £? = a v B t cx -(SI ■ A^)Bi (see 
where A(") determines the inhomogeneity of turbulence. 
The S term in the electromotive force caused by the 
shear-current effect is given by £f = -((5x(VxB)), cx 
— (W • V)Bi (see H3), where the 5 term is proportional 
to the mean vorticity W = VxU which is caused by the 
differential rotation. 

During the generation of the mean magnetic field in 
both cases (in the aSl dynamo and in the shear-current 
dynamo), the mean electric current along the original 
mean magnetic field arises. The a effect is related to 
the hydrodynamic helicity cx (fl ■ A^) in an inhomoge- 
neous turbulence. The deformations of the magnetic field 
lines are caused by upward and downward rotating tur- 
bulent eddies in the aft dynamo. Since the turbulence is 
inhomogeneous (which breaks a symmetry between the 
upward and downward eddies), their total effect on the 
mean magnetic field does not vanish and it creates the 
mean electric current along the original mean magnetic 
field (see 0). 

In a turbulent flow with the mean differential rotation, 
the inhomogeneity of the original mean magnetic field 
breaks a symmetry between the influence of upward and 
downward turbulent eddies on the mean magnetic field. 
The deformations of the magnetic field lines in the shear- 
current dynamo are caused by upward and downward 
turbulent eddies which result in the mean electric current 
along the mean mag netic field and produce the magnetic 
dynamo (see [25l |26| ) . 

Note that the differential rotation is described by 
the gradient tensor of the mean velocity field ViUj = 
(dU)ij + £y„(VxW)„/2, where the symmetric part of 
the gradient tensor (dU)ij — (V^t/j + \7jUi)/2 is given 
by 

(9U)ij = ^[(e.xr^V, + (e z xr) 3 V,] (60.) , (23) 

and the mean vorticity W in cylindrical coordinates is 
given by 

W = -p(e p W z -e z W p )(Sn) , (24) 

and in spherical coordinates the mean vorticity is 

W = rsin0(e r V0-e e V r )(5fi) . (25) 

The nonlinear coefficients <7q(B) and <J\(B) defining 
the shear-current effect and the nonhelical a effect are 
determined by Eqs. (|A64|) and l|A65|) in Appendix A. 
The nonlinear dependencies of the parameters cro (B) and 
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FIG. 1: The nonlinear coefficient ao(B) defining the shear- 
current effect for e = (thin solid) and for e = 1 (thin dashed) ; 
and the nonlinear coefficient o\{B) defining the nonhelical a 
effect for different values of the parameter e: e = (thick 
solid); e = 17/21 (thick dashed-dotted); e — 1 (thick dashed). 

Ci{B) are shown in FIG. 1 for different values of the pa- 
rameter e. The background magnetic fluctuations caused 
by the small-scale dynamo and described by the param- 
eter e, increase the parameter ao(B). For a weak mean 
magnetic field B <C B eq /4 the parameter ao(B) is given 
by a (B) = (4/45) (2 - q + 3e) (see HH), where q is the 
exponent of the energy spectrum of the background tur- 
bulence. The latter equation is in agreement with that 
obtained in |25j| where the case a weak mean magnetic 
field and e = was considered. In this equation we ne- 
glected small contribution ~ 0[(4B / B cq ) 2 }. The mean 
magnetic field is generated due to the shear-current ef- 
fect, when <7q(B) >, i.e., when the exponent of the energy 
spectrum q < 2 + 3e. Note that the parameter q varies 
in the range 1 < q < 3. Therefore, when the level of the 
background magnetic fluctuations caused by the small- 
scale dynamo is larger than 1/3 of the kinetic energy of 
the velocity fluctuations, the mean magnetic field can be 
generated due to the shear-current effect for an arbitrary 
exponent q of the energy spectrum of the velocity fluctu- 
ations (see H^l). 

For the Kolmogorov turbulence, i.e., when the expo- 
nent of the energy spectrum of the background turbu- 
lence q — 5/3, the parameters ao(B) and ci(-B) for 
B < B_ cq /4 are given by a (B) = (4/135) (1 + 9e) 
and cn{B) = (2/135) (17- 21e). For B > B cq /4_they 
are given by a {B) = -(11/135) (1 + e) and a x {B) = 
(2/135) (1 + e). It is seen from these equations and from 
FIG. 1 that the nonlinear coefficient (Jq{B) changes its 
sign at some value of the mean magnetic field B = B*. 
For instance, B* = 0.6B eq for e = 0, and B. t = 0.3B C q 
for e = 1. However, there is no quenching of this ef- 
fect contrary to the quenching of the nonlinear a effect, 
the nonlinear turbulent magnetic diffusion, the nonlinear 
Ox J effect, etc. 



FIG. 2: The nonlinear coefficient 8q(B) defining the f2xj 
effect for different values of the parameter e: e — (dashed) ; 
e = 2/7 (dashed-dotted); e = 1 (solid). 



The mean differential rotation causes the nonhelical a 
effect, W*ai(B)V z (5n) [see Eqs. and (Q3J)], which 
is independent of a hydrodynamic hclicity. It follows 
from the asymptotic formula for cr\{B) at B ^> £? oq /4 
that there is no quenching of this effect contrary to the 
quenching of the regular nonlinear a effect (see Section 
III- A). These two kinds of the a effect have opposite 
signs. Thus, the total a effect should change its sign 
during the nonlinear growth of the mean magnetic field. 
The nonhelical a effect vanishes if the mean rotation is 
constant on the cylinders which are parallel to the rota- 
tion axis. Note that <J\{B = 0.1£? cq ) = for e = 1. 

The 5 term in the electromotive force which is re- 
sponsible for the shear-current effect has been also cal- 
culated in pfl l4l| for a kinematic problem using the 
second-order correlation approximation (SOCA). How- 
ever, these studies did not found the dynamo action in 
nonrotating and nonhelical shear flows. Note that the 
second order correlation approximation (SOCA) is valid 
for small hydrodynamic Reynolds numbers. Indeed, even 
in a highly conductivity limit (large magnetic Reynolds 
numbers) SOCA can be valid only for small Strouhal 
numbers, while for large hydrodynamic Reynolds num- 
bers (fully developed turbulence) the Strouhal number 
is unity. Our recent studies for small hydrodynamic and 
magnetic Reynolds numbers (using spectral r approxima- 
tion) also did not found the dynamo action in nonrotating 
and nonhelical shear flows in agreement with |40l l4l| . 



C. The nonlinear coefficient S$(B) denning the 
f2xJ effect 

The 8 term in the electromotive force which is caused 
by a uniform rotation, describes the fl x J effect. This 
effect in combination with the differential rotation can 
cause a generation of the mean magnetic field even in a 
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nonhelical turbulent flow (see |42j |43|, |4J, |45j and |2Jj), 
where J is the mean electric current. The nonlinear co- 
efficient Sq(B) defining the fix J effect is determined by 



* 4 {Ci + C 3 } y - (I - e)(* a + 4* 3 ){^i 
+A 2 } y + (l + e)(A^\y)+A^(y))\ 



y=iB 



(26) 

where the functions ^k{X}y are determined by 
Eqs. (|A44J| in Appendix A. The parameter 5q(B) is de- 
termined by the contributions from the (5(B) term, the 
r)ij(B) term and the Ky/-(B) term in the nonlinear elec- 
tromotive force caused by a uniform rotation. The 
nonlinear coefficient 6q(B) is shown in FIG. 2 for differ- 
ent values of the parameter e. The asymptotic formulas 
for the coefficient 5ff(B) for a weak mean magnetic field 
B < ,B cq /4 are 



6$(B) 

and for £? ^> B cq /A are 



135 



(2-7e) 



3^(34+19.) 



(27) 



(28) 



Asymptotic formula l|27|l for a weak mean magnetic field 
(£? <C .Beq/4) coincide with that obtained in [2j| for g = 
5/3. 



D. The dynamical equation for the function x c (B) 

The function x c (B) entering the magnetic part of the 
a effect [see Eq. ()15|l] is determined by the dynamical 
equation 



[£■6 + v • :f (x) 1 



dt \i 0/ 

-V • (U X C ) - X c /T 



(29) 



(see, e.g., 0,1^), where is the nonadvective flux 

of the magnetic helicity which serves as an additional 
nonlinear source in the equation for x c ( see HH 122^1. 
Ux c is the advective flux of the magnetic helicity, U 
is the differential rotation, and T = (l/3)(lo/L) 2 Km is 
the characteristic time of relaxation of magnetic helic- 
ity. Equation l|29|) was obtained using arguments based 
on the magnetic helicity conservation law. The func- 
tion x° is proportional to the magnetic helicity, x° = 
2x m /(9^77 T p ) (see QjJ), where x™ = ( a • b) is the mag- 
netic helicity and a is the vector potential of small-scale 
magnetic field. The physical meaning of Eq. (|29|l is that 
the total magnetic helicity is a conserved quantity and if 
the large-scale magnetic helicity grows with mean mag- 
netic field, the evolution of the small-scale helicity should 



somehow compensate this growth. Compensation mech- 
anisms include dissipation and various kinds of transport 
(see [Uli^). 

In order to demonstrate an important role of the non- 
advective flux of the magnetic helicity, let us consider 
a local model in cylindrical coordinates, when the mean 
magnetic field depend only on the vertical coordinate z 
and A' > A/r, where A' = dA/dz and B = Be v - A'e r . 
Since 



dA 



dB 
~dt 



= EL- DA' 



(30) 
(31) 



we obtain that 



Then in steady state Eqs. and (J22J) yield 
Z 2 r - 2DT [ z x) = const , 



(32) 



(33) 



where £ r = rj B B'. Here we neglected the last term in 
Eq. H29|) which, e.g., for galactic dynamo is very small. 
In a steady-state for fields of even parity with respect 
to the disc plane, we obtain the solution of Eq. 1)33(1 for 
positive C D 



dB = v /2|C r D| / V\rW\d~z , (34) 
\F{B)\ J W 



where J r i x) = C\J 7 (B)\\ X V (z)\. The crucial point for the 
dynamo saturation is a nonzero flux of magnetic helicity. 
It follows from Eq. I)34|) that this saturation mechanism 
is nearly independent of the form of the flux of magnetic 
helicity. In that sense this is a universal mechanism which 
limits growth of the mean magnetic field. If we assume 
that 1^(5) I ~ _B -27 *, we obtain that the saturated mean 
magnetic field is 



\CD\ 



V\x v (z)\ dz 



B,, 



(35) 



where we redefined the constant C, we took into account 
that w B (B) cx B cq /B for B ^> B cq /4, and we restored 
the dimensional factor B cq . Note that the nonadvective 
flux of the magnetic helicity was chosen in in the 

form T (x) = Cx v 0v{B)B 2 r] ( J;\B)(Vpo)/po. This corre- 
sponds to 7* = I in the function For the specific 
choice of the profile — sin 2 (7T2i/2) we obtain 



— ^\CD\B cq co S ( — 



1 + e - f 7TZ 

— ; -B rn tan — 

4|iL| cq V 2 



(36) 
(37) 
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where we have now restored the dimensional factor B eq . 
The boundary conditions for B v are B v {z = 1) = 0, 
B'Jz = Q) = 0, and for B r are B r {z = 1) = 0, 
B r (z = 0) = 0. Note, however, that this asymptotic anal- 
ysis performed for B 3> B cq /4 is not valid in the vicinity 
of the point z = 1 because B(z = 1) = 0. 



E. The dynamo waves 

In order to elucidate the new effects caused by the dif- 
ferential rotation, let us consider first a kinematic prob- 
lem in a spherical geometry. Following we study 
dynamo action in a thin convective shell, average the lin- 
earized equations l|13|l and (|14|) for A and B over the 
depth of the convective shell. Then we neglect the curva- 
ture of the convective shell and replace it by a flat slab. 
These equations are obviously oversimplified. However, 
they can be used to reproduce basic qualitative features 
of solar and stellar activity (see, e.g., ^3|)- We are in- 
terested in dynamo waves propagating from middle solar 
latitudes towards the equator. We seek for a solution of 
the obtained equations in the form of the growing waves, 
A, B oc exp(jt) exp[i(u)t — K-R)], where the growth rate 
of the dynamo waves with the frequency 



-ai 



D\S K \ sgn(S K ) 



(38) 



is given by 



7 : 



'D\S K \ 



07 



K z 



(39) 



The frequency and the growth rate of the dynamo waves 
are written in a dimensionless form. Here 



r>* # k z 



ai = a + W* o\ S z , 

(Tl = W* (To Sk 

s k = K e v r {sn) - K r v 9 (sn) , 

S z = cos 9 V r (5fl) -sin0V e (tfft) , 
K z = cos 9 K r — sin 9 Kg . 

The total a effect, a;, is a sum of the usual a effect 
(caused by helical motions) and a nonhclical contribu- 
tion, W* ai S z , due to the effect of the the mean differen- 
tial rotation on the small-scale turbulence. The parame- 
ter ai describes both, the shear-current effect determined 
by W* do Sk term, and the fix 3 effect determined by 
f2* 5q K z term. Even in nonhelical turbulent motions, 
the mean magnetic field is generated due to the shear- 
current effect and the Ox J effect. 



IV. DISCUSSION 

Let us discuss the nonlinear effects. It was shown re- 
cently in [2(j that the algebraic nonlinearity alone (i.e., 



algebraic quenching of both, the a effect and turbu- 
lent magnetic diffusion) cannot saturate the growth of 
the mean magnetic field. Note that the saturation of 
the growth of the mean magnetic field in the case with 
only an algebraic nonlinearity present can be achieved 
when the derivative of the nonlinear dynamo number 
with respect to the mean magnetic field is negative, i.e., 
dD N {B)/dB < 0. HereL>Ar(i?) = a(B)/[r] A (B)r] B (B)} is 
the nonlinear dynamo number. Thus, when the nonlinear 
dynamo number decreases with the growth of the mean 
magnetic field, the nonlinear saturation of the magnetic 
field is possible. 

In this study we showed that the differential rotation 
of fluid can decrease the total a effect. In particular, the 
mean differential rotation causes the nonhelical a effect, 
W* 0i (B) V z (5f2), which is independent of a hydrody- 
namic helicity. We demonstrated that there is no quench- 
ing of this effect contrary to the quenching of the regular 
nonlinear a effect, a(B) = x v ^ v (B) + a n + x c (B)(j) m (B). 
In this study we found that these two kinds of the a ef- 
fect have opposite signs. Thus, the total a effect should 
change its sign during the nonlinear evolution of the mean 
magnetic field, and there is a range of magnitudes of the 
mean magnetic field, where the nonlinear dynamo num- 
ber decreases with the growth of the mean magnetic field. 
Therefore, the algebraic nonlinearity alone can saturate 
the growth of the mean magnetic field if one take into ac- 
count the effect of differential rotation on the nonlinear 
electromotive force. For instance, the nonhclical a effect 
causes a saturation of the growth of the mean magnetic 
field at the base of the convective zone at B < 2B cq 
(see below), where B eq is the equipartition mean mag- 
netic field. However, the nonhclical a effect vanishes if 
the mean rotation is constant on the cylinders which are 
parallel to the rotation axis. 

In this study we also demonstrated that the mean dif- 
ferential rotation which causes the shear-current effect, 
increases a growth rate of the large-scale dynamo insta- 
bility at weak mean magnetic fields, and causes a sat- 
uration of the growth of the mean magnetic field for a 
stronger field. 

The nonlinear shear-current effect and the nonhelical 
a effect become very important at the base of the con- 
vective zone (see below). When we apply the obtained 
results to the solar convective zone, we have to take into 
account that all physical ingredients of the dynamo model 
vary strongly with the depth H below the solar surface 
and we have to use some average quantities in the dy- 
namo equations. We use mainly estimates of govern- 
ing parameters taken from models of the solar convec- 
tive zone (see, e.g., I n particular, in the upper 
part of the convective zone, say at depth H ~ 2 x 10 7 cm, 
the magnetic Reynolds number Rm ~ 10 , the maximum 
scale of turbulent motions lo ~ 2.6 x 10 7 cm, the char- 
acteristic turbulent velocity in the maximum scale Iq of 
turbulent motions uq ~ 9.4 x 10 4 cm s _1 , the fluid density 
po ~ 4.5 x 10~ 7 g cm" 3 , the turbulent magnetic diffu- 
sion r] T ~ 0.8 x 10 12 cm 2 s _1 and the equipartition mean 
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magnetic field is B cq = 220 G. Thus, in the upper part 
of the convective zone the parameters W* ~ 10~ 3 — 10~ 4 
and f2* ~ 5 x (10~ 3 — 10~ 4 ). According to various mod- 
els, the ranges of the dynamo number D « 10 3 — 10 6 can 
be considered as realistic for the solar case. At the base 
of the convective zone (at depth H ~ 2 x 10 10 cm), the 
magnetic Reynolds number Rm = IqUo/t) ~ 2 • 10 9 , the 
maximum scale of turbulent motions Iq ~ 8 x 10 9 cm, the 
characteristic turbulent velocity uq ~ 2 x 10 3 cm s _1 , the 
fluid density p ~2x 1CP 1 g cm~ 3 , the turbulent mag- 
netic diffusion r\ T ~ 5.3 x 10 12 cm 2 s -1 . The equipartition 
mean magnetic field B cq = 3000 G. Thus, at the base of 
the convective zone the parameters W* ~ 1 — 10 and 
f2* ~ 5 — 50. Thus, the effects of the differential rotation 
(the nonlinear shear-current effect and the nonhelical a 
effect) become very important at the base of the con- 
vective zone. Since these effects are not quenched, they 
might be the only surviving effects. 



M 



ijpq 



(l/2k' 2 )(k rn S n p + k n 5< 



— Oh h 
mp a™ ^mnp J i 







where Dij(k) = 2sij m k m (k ■ f2)/fc 2 . The tensors 
^mn(U), ££ mn (U) and Jg m „(U) are given by 



J" 7 . (U) = 



$im$jq$np 



E" 



ijrnn 



t (U) 
(U) 



.9 1 

^iq^jn^rnp ~r~ ^im^jn^q ^pUq , 

^^jq^im^pri ^im^pn^jq "t - ^jn^pm^iq 
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APPENDIX A: EFFECTS OF UNIFORM AND 
DIFFERENTIAL ROTATIONS 

The method of the derivation of equation for the non- 
linear electromotive force in a rotating turbulence is sim- 
ilar to that used in [2(| for a nonrotating turbulence 
with an imposed mean velocity shear. In the frame- 
work of a mean-field approach we derive equations for 
the following correlation functions: /y(k) = L(Ui,Uj), 
/iij(k) = L{bi]bj) and <ftj(k) = L(bi;v,j), where L(af,Cj) 
is determined by Eq. (JSJ. In order to exclude the pres- 
sure term from the equation of motion Q3J we calculate 
Vx(Vxu). Then we rewrite the obtained equation and 
Eq. |0J in a Fourier space. The equations for these cor- 
relation functions are given by 



at 



dh tj (k) 
dt 



%j( k ) 
dt 



+i(k-B)$. 



tjrnn \ 
(M) , 



Ft 



t(U)A m „-i(k.B)$ 



i jinn 

h < bh N 



-i 



i j 



(Al) 



(A2) 



-/ly'(k) 



' tjmn 



I 9 . 



Dg 



N 



where the mean 
ferential rotation, 



velocity U describes the dif- 
= 0«(k) - Sfj»(-k), 



Fy(k) = (F(k) % -(-k)) 



i(k)F.,(-k)>, F(k) = 



kx (kxF(k))/fc 2 /9 . The tensors M-j q and are given 
by 



M 



12 

ijpq 



D ip (\Li)6jq + D jq (\s. 2 )Sip 

Ai(k 2 ) 



M + Af Si 

tjpq ' %3pq 



1)9. 



t 2^1 rn k rnn (^£ip n 5jq -\- Siqn^ip) 



(see [23, Hg), where is the Kronecker tensor, fey = 
kikj/k 2 . Equation (|A1|) - (|A3|) are written in a frame 
moving with a local velocity U. For the derivation of 
Eqs. l|A"T ]) -l|A3 j) we used the relation 

£-ijnQnk {^inlkj £jnlki)kn^l — ^ijn^n(k ' ^) j 

which applies to arbitrary vectors k and CI (see psj - ). 
The source terms if- , A and I?- (which contain the 
large-scale spatial derivatives of the mean magnetic field 
and the second moments) are given by 



iL = 



m = 



1?. 



^(B-V)^P + [g qj (k)(2P m (k)-6 m ) 



+g qi (-k)(2P jn (k) - 5 jn )}B n>q 
1 



B k <i> (P ^ 



(A4) 



(B-V)^P - [g iq (k)6 3n + g jq (-k)S in ]B niq 



-R k $>( pS> 
O n ,q^n^ij q 



(A5) 



-(B-VX/ij + hij) + h iq {2P jn {k) ~ 5 3n )B n , q 



ijq -r hij q ) 
(P) 



(A6) 



(see 123), where V = d/dK, ^'(k) = gij (k)+ 9ji (-k), 
and Bij = Vji?i, fff, hfj and gfj are the third mo- 
ments appearing due to the nonlinear terms, fij q — 

(A3) (1/2)8 fij /dk q , and similarly for hij q and $\f q • To derive 



Eqs. {HJl-JIlJl we used the identity: 



dK dQ (k p + K p /2)Bp(Q) exp(iK-R) 



x (m(k + K/2 - Q)u 3 -(-k + K/2)) 



i(k-B) + -(B. V) 



f l3 (k^)-\kp d -M^B 



(A7) 



(see We took into account that in Eq. (|A3|> the 

terms with symmetric tensors with respect to the indexes 
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"i" and "j" do not contribute to the electromotive force 
because £ m = e m jigij- In Eqs. HA1 M A3 |I we neglected 
the second and higher derivatives over R. To derive Eqs. 
(|A1|) - (|A3|) we also used the following identity 



iki I fij(k - -Q, K - Q)U P (Q) exp(iK-R) dK dQ 



i ( dfij 



4 V dk s 



(A8) 



(see [23). We split the tensor of magnetic fluctuations 



into nonhelical. 



hij, and helical, 



parts. The helical 



part of the tensor of magnetic fluctuations depends on the 
magnetic helicity and it is not determined by Eq. i|A2|l . 

The tensor is determined by the dynamic equation 
which follows from the magnetic heli city conservation ar- 
guments H H3 (see also [H Q E El IH HI ) . 

First, we consider a nonrotating and shear free turbu- 
lence (SI = 0; V.;U = 0), and we omit tensors I£ TOn (tJ), 



E: 



1 j m il 



(U) and J[ jmn (V) in Eqs. 



l|Al )l -l|Al )l . First we 

solve Eqs. ijAl(l - HA3ll neglecting the sources 7^-,/^-,J?- 
with the large-scale spatial derivatives. Then we will 
take into account the terms with the large-scale spatial 
derivatives by perturbations. We start with Eqs. (|A1|) - 
(|A3|) written for nonhelical parts of the tensors, and then 
consider Eqs. HA1[) - I|A3[I for helical parts of the tensors. 

We subtract Eqs. (|A1I) - (|A3|) written for background 
turbulence (for B = 0) from those for B ^ 0, use the 
t approximation [which is determined by Eqs. ©-(JSJ 
neglect the terms with the large-scale spatial derivatives, 
assume that rjk 2 <C r~ 1 (fc) and vk 2 <C r^ 1 (k) for the 
inertial range of turbulent fluid flow, and assume that 
the characteristic time of variation of the mean magnetic 
field B is substantially larger than the correlation time 
r(k) for all turbulence scales. We split all correlation 
functions into symmetric and antisymmetric parts with 



respect to the wave number k, 



fij — f. 



0) 



where f\ 
and/, 



e -9-, J13 J%3 ■ jij 



f, 



(a) 



(a) 



= [fij(k) + fij(— k)]/2 is the symmetric part 
[/ij(k) — fij(— k)]/2 is the antisymmetric part, 



and similarly for other tensors. Thus, we obtain 



/?(k) 



Y^[(i + Wff' s) (k) + ^ M ( k )]' 



'(k) 



T ^-[^ s) (k) + (i + ^ ) (k)], 

(A10) 

ir(k-B) rj(0|a)ri ^ (0i .) 



(All) 



neglect the large-scale spatial derivatives, i.e., they are 
proportional to the first-order spatial derivatives. 

Now we take into account the large-scale spatial deriva- 
tives in Eqs. (| Al|> - (| A3|> by perturbations. Their effect 
determines the following steady-state equations for the 
second moments fij, hij and <?jj-: 



$°(k) 



^■ } (k) 



/^ ) (k)+^r(k.B)$ ^ ( f s) (k)+r4 



^(k)- i r(k.B)$<f s) (k)+r4, 



(A12) 
(A13) 



(k.B)(^ a) (k)-^ ) (k))+r/?. , 



(a), 
ij 



(A14) 



where $ 



(M,s) 



^r ) (k)+c ) (- 



■k)]/2. Here L 



and cjij denote the contributions to the second moments 
caused by the large-scale spatial derivatives. The corre- 
lation functions of the background turbulence /^'"•'(k) 
and hf*' a \\t) are determined by the inhomogeneity of 



turbulence [see Eqs. 
(|AT2 | - ljAT4| yield 



© and (HOI)]. The solution of Eqs. 



4> 



Mf,) f u ^r(k-B) 
W ~ 1 + 2^ 



(f(0,a) _ J,((0,a))\ 
\J ran ran } 



x(l + 2?/?)((5 n j(5 m fc 5 m jS n k + k n ^8 m j fcmfc^nj) 

T B jtk 



-2(e+ 2lp)(k n jS m k - kmjSnk) 



], The correlation functions f-j (k), h\j'(k) and g\j'(\s) are 
of the order of ~ 0(V 2 ), i.e., they are proportional to 
the second-order spatial derivatives. Thus fij + fij is the 
nonhelical part of the correlation function of the velocity 
field for a nonrotating turbulence, and similarly for other 
second moments. 

Next, we solve Eqs. l|Al|l - l|A3|l for helical parts of the 
tensors for a nonrotating turbulence u sing the same ap- 
proach which we used before (see also |26|). The steady- 
state solution of Eqs. ijAlfl and l|A3(l for the helical parts 
of the tensor reads: 



(1 + 2V-) 2 



+ 

(A15) 



(M,H) 



(k) 



2<r(k.B) (0lH) 



where $^' H ^(k) 



(H) 

: 9ij 



1 + V 
(k) 



(A16) 



<4 H) (-k) and/f H) (k) is 



(A9) the helical part of the tensor for velocity field of the back- 
ground turbulence. The tensor h\j is determined by the 

dynamic equation Since f^ and h\f are of 

the order of O(V) we do not need to take into account 
the source terms with the large-scale spatial derivatives. 
Now we determine the nonlinear electromotive force 



£i(r = 0) = (l/2) £mm /[$^ ff) (k) + $^(k)]dk in a 



(see |2fJ), where f tj , h l3 and fa are solutions without the nonrotating and shear free turbulence: 



sources I(j , Ifi 



and Ifj, 



V(k) = 2(rk-B) 2 . The correla- 



tion functions f^ J (k), h\^(k) and g^' (k) vanish if we 



ir(k-B) (0 H) _ (H) 

-1 1 / \J ran ran I 

l + V 



1 + 2^ 



{l, 



9 

ran 
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+i(k-B)[/, 



(0,o) 
mn 



u (0,o) , irf _jh \|i 
mn 1 rmn mn/JJ 



(A17) 



To integrate in k-space in the nonlinear electromotive 
force we specify a model for the background turbulence 
[see Eqs. (|9"j)- (ll0[l ]. After the integration in k-space we 
obtain the nonlinear electromotive force: 



(A18) 



dBi/dRj, Eijk is the Levi-Civita tensor, 



where B 

and the tensors a^- and byj. are given by 



1 

6 T ° 



A^ ( y/20) E ijn - 4 X) (V2/?)e mm A 



xV„[(u 2 )(°) - (b 2 }(°) 
+X c (B)0" l (/3)]% , 



4>l{B) Eijk + 4>2{B) Eij n Pnk 



~4>?,{B) Sink Pnj 



(A19) 



(A20) 



where ftj = BiBj/B 2 , the quenching functions 4> v (f3), 
<j? n {[i) and fc (B) are determined by Eqs. JT7J) , CHJ) 
and respectively, /3 = 4B/ (u y / 2[ip), e = 

(b 2 )(°) /(u 2 )(°\ and all calculations are made for q = 5/3. 
The parameter \ v — — T o ft v /3 is related to the hydrody- 
namic helicity fi v of the background turbulence, and the 
function x c (B) = (r/3/xp)(b ■ (Vxb)) is related to the 
current helicity. These parameters are written in the di- 
mensional form. To integrate over the angles in k-space 
we used the following identity: 



1 + a cos 2 9 



d9dip = Arfy + A 2 (3 l3 , (A21) 



where a = f3 2 /r(k), and 



2tt 
a 



(a + 1) 



arctan( v / a) 



, 2tt 
A 2 = 



(a + 3) arctan(^)_ 3 



(A22) 
(A23) 



(for details, see [H IH). The functions a1 1} (/3) are given 
by 

A£H0) = ^£^f±dX, (A24) 

where A 2 = /? 2 (fc/fco) 2 / 3 = a, and we took into account 
that the inertial range of the turbulence exists in the 
scales: Id < r < lo- Here the maximum scale of the tur- 
bulence lo <C L, and /<; — lo/Re 3 ^ 4 is the viscous scale 
of turbulence, Re = IqUq/v is the Reynolds number, v is 
the kinematic viscosity and L is the characteristic scale 
of variations of the nonuniform mean magnetic field. For 



very large Reynolds numbers kd — I J 1 is very large and 
the turbulent hydrodynamic and magnetic energies are 
very small in the viscous dissipative range of the tur- 
bulence < r < Id- Thus we integrated in A n over k 
from fco = Iq 1 to oo. We also used the following identity 
/„ A n (a(?))?d? = (2tt/3)a1 1) (/3). The explicit form of 

the functions Ak((3 2 ) and A^\(3), and their asymptotic 
formulas are given in j2(|. 



We use an identity Bj^ = 
which allows us to rewrite Eq. 
tive force in the form of Eq. J5J 



(dB) tJ + £y„(VxB)„/2 



.181) for the electromo- 
whcre 



/(B) 



1/ 
1 



+ aji ), Vf(B 



1 

2 £ kji a ij 



TJij (B) ^(Eikpbjkp ~1~ E-jkpbikp) 



(A25) 
(A26) 



8i = 



), 



jfc(B) 



- --^{bijk + hkj). 

(A27) 

(see [13). Using Eqs. (|X25f - jATTf and PT9f - jJ20)) we 

derive equations for the coefficients defining nonlinear 
electromotive force for a nonrotating turbulence. In par- 
ticular, 



V 



err 



[X w <ra+X C (B)<T (/?)]% , (A28) 



V d (B 

2(B) 



where 



Vij 



V d (B) 



B 2 

rj T <pi(B)6ij , 



(B • V)B 



(A29) 
(A30) 



1 



MB)(A^ -eAW) , (A31) 



and A< g ) = (V B 2 )/B 2 . Note that Eqs. EH2l) - (IX20l) 
and HA28(I - (|A30() for a homogeneous and nonhelical back- 
ground turbulence coincide with those derived in |26| . 

Now we study the effect of a mean uniform rotation 
of the fluid on the nonlinear electromotive force in a 
shear free turbulence. We consider a slow rotation rate 
(rfi <C 1), i.e., we neglect terms ~ 0(Q 2 ). We also neglect 
terms ~ 0(V 2 ). However, we take into account terms 
~ 0(f2jVj), that is possible by the following symmetry 
reasons. The tensor r^Vj is a pseudo tensor, while f^ll, 
and V;Vj are true tensors. This implies that a pseudo 
tensor quantity includes terms oc Cl^Vj, but does not in- 
clude terms oc f2jfij and oc VjVj. On the other hand, 
a true tensor quantity does not include terms oc Q,iVj, 
but it may include the terms oc f2jf2j and oc ViVj. 
The steady-state solution of Eqs. (|A1|) and (|A3I) for the 
nonhelical parts of the tensors for a rotating turbulence 
reads: 



-ij 



-{z(k-B)^ 

-{»(k.B)[/«(k) 



if.} 



(A32) 



(k)] 



(A33) 
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where N{ jpq (fl) = 5 lp 8 3q - rM% pq and iV£(fi) = 5 {j - 
tD%. Here we use the following notations: the total cor- 
relation function is fy = fij + f°j, where /y = + /y 
is the correlation functions for a nonrotating turbulence, 
and ffj determines the contribution to the correlation 
function of the velocity field caused by a uniform ro- 
tation. The similar notations are for other correlation 
functions. Now we solve Eqs. lA2jl . 1A32|) and (|A33j) by 
iteration which yields 



L ijpqJ PQ 



(A34) 



-r{i{k-B)<\^' M /£"($)}. (A35) 



7g(k) - r{^ B ft„+<(k.B)[/n_^g] 

+^ n (/M)}> 

(k) = ffy(k) - gfii-k), the source terms 



(A36) 



^'"(sg) and ^'"(/S^g) are determined by 
Eqs. (|A4J) - ljA6J) . where /y, /ly, are replaced by 
hfj, gV, respectively. The solution of Eqs. HA34jl - ljA36|) 

yield equation for the symmetric part $>j w,f2,s -' of the 
tensor: 



(m,o,. 



'(k) = 



1 + 20 



{D%g in - D? n g jn + ir(k-B)[lff 



ijpqjpql 



(A37) 



Thus, the effect of a uniform rotation on the 
nonlinear electromotive force, £p(r = 0) = 
(l/2)e inm J $mn n ' s ^ dk, is determined by 

Si = ernmj {D%g mp + ir(k-B) [M% npq f pq 

+If£-I^]+I%£}dk. (A38) 
Now we use the following identities: 

£-inniD np (ji n p — < 2^l rn (^ki m g pp k nm g n ^ , 

^■inra^ np9mp — 2£ (T'igpp ^ndni) , 
£inmM mn pqf p q — ^^m^nmifin fni) ; 
i£inmM mn pqf p q = A(Tif pp T n fi n ) , 



where Tj = r2 m T m i p V p . We also take into account that 



knfni — /2)V 'nfni 7 ^n/m /2)V nfin : 

k n g n i = (*/2)V n g n i . 

These equations follow from the condition Vu = 0. 
Thus we obtain that the effect of a uniform rotation 
on the nonlinear electromotive force is determined by 
Si=^B j + b% k B j<k , where 



2t O rn 
1 + 2-0 

l - 20 



[(u 2 )(°) - (b 2 }(°)] 



2(ir 



2\(0) 



1 + 20 

+ k im S jk [(u 2 >(°) + (b 2 )(°)]Uk 



(A39) 



,n _ 

u ijk — 



2t Qrn 

(1 + 20>) 2 

20 fcjm ^ifc + 



2(V-1) + 



1 + 2-0 



1 + 20 



kijmk 

hj8 mk [(u 2 )(°) 
d 



-(b 2 )(°)]-(l + 2^)fc% m ^-[(u 2 )W 



-(b 2 )(°)]Uk 



(A40) 



and we used the identities: 



(k-B)V„0 = 2-ipkjB.j 



k; 



|g = 4r 2 (k.B)S t -2 ((Z -l)^ 
Now we use the following identities: 



BjKij mn Q m A n = {(Ci + C 2 ) [OjAj + fijAj 

+ 5 y (n-A)] + [Ci + 3C 3 ) <5 y - (n-j§)(A-3)}Bj , 
KijQj(B-A) = [A, CkAj + A 2 6 Z] (ri-/3)(A-/3p, , 
B 3 K l3 = (A 1 +A 2 )B i , 

B 3 ,i K jm n m = A 1 v,(rj-B) + ~ a 2 Ai B) (n.B) , 
% m i^n m = [A! (n- v) + \ A 2 (fi.A( B >)] ^ , 



Bj, k Kij mk Q, m — Cj (f2j <5jfc + 0^ (5.^; j-^.n, . ^ 



) k + I {c 3 [2Aj B) (n.B) + ^(A^'-B)] + 



(ftA( B )) 



+ --C3.v)(n.B) 

5j (B* V) Kijjjin^lm 



C 2 (0./3)(A( B )./3)B 4 | - C 3 i (0-/3) (Bx(VxB)) , 



(Cx + C 3 ) {[A\ B) (n-B) + ^(A^-B) + B,(n-A (B »)]B 2 - 2[%B 2 
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+ <5 ij (0-B)](Bx(VxB)) j } + {C 2 + 3C 3 ) (n-B)(A^-B) , 

where 

/k ' ' sin — — 

l _|_ a cos 2 Q d8 dip — Ci(SijS mn + 6i m Sj n + Si n Sj m ) + C 2 fiijmn + C 3 {8ij(3 mn + 5i m (3j n 

~\~3inPjm ~t~ fijmPin "4~ ^jnPim ^mnfiij) , (-^-41) 



and 



0l ~ 2V 



^ _l_ -^2 arctan(\/a) 5a 



1 

3 



I 

C 2 = A 2 - 7A 1 + 35C*i , C 3 = A 1 - 5(7i 
Integration in k-space yields 



°ijk 



= ^2| jB30 . A (B) + ^ 4Q . A (B) + ^_ ^ (^ 5 (n.A( B ))B 2 + [£ 6 (B.A( B >)+ £V(B-V)](n-B) 
+£ 8 n-(Bx(VxB))j -£ 9 £ iim (VxB) m (n-B) + E w %a| u) + E n fi^" ] + % ^ 12 (n-A<">) 

+E 13 (n-B)(B-AW)j + e £ 14 fy-Af > + £7 1B n^f + % (# 16 (fl-A^) + E 17 i(n-B)(B-A( 6 ))j | , 

(A42) 
(A43) 

I 



2 



where 

^{X}, = 3XW(y)-|-X(y 2 ), 

* 2 {X}, = 4X( 2 )(y)- J-X(y 2 ), 
* 3 {X} y = 6X( 2 )(y)--X( 2/ 2 ) + 

* 4 {X} y = [2-(l + e)(2 g -l)]* 2 {X}„ 

+4(l-e)* 3 {X} a -(l + 3e)X( 2 )(y) , 
* 5 {X}„ = 2* 2 {X}, + X( 2 )(y) , 
^ 6 {X} y = -2M> 2 {X} y + X^(y) , (A44) 

X' = dX/dz, and all calculations are made for q = 5/3, 



Si 
E 2 

E 3 



* 4 {Ci}„ + (1 - e)* 2 {Ai} w + 2eAf\y) 
* 4 {C 1 },-(l-e)(vI/ 2 + 4* 3 ){ J 4 1 }, 
+ (l + e)4 2) (y) 



Ev — —Ez -\- — 
7 2 2 



y=V2/3 

(1 - e)* 2 {A 2 }„ 



+2e4 2) (y) 



y=v / 2/9 



£4 
#5 

£-6 
#7 

E$ 

Eio 
E\ 2 



E4 + - 



(l + e )4 2) (y)-(l- £ )(* 2 



+4* 3 ){A 2 } y 



v=y/20 

* 4 {C 2 } y - 2(1 + e) y 2 *i{C 2 + 3C 3 h 
Eg — Ei — Es , 
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£7 14 

-^16 
£17 



*4{C 3 } y=V2/3 

4(l + e)/3 2 vI/ 1 {C 1 +C 3 }^ /3 , 

* 5 {Ci + C 3 }^ , E n = E w + A?\V2(3) , 
Ew-^siAi+A^^ , 

^f 5 {C 2 + 3C 3 } y + A^(y) 

ly=V2P 

* 6 {Ci + C 3 }^ fj , £ 15 = E 14 + A {2) (V2/3) , 
E 14 - ty 6 {A 1 +A 2 } V 5p , 

^ 6 {C 2 + 3C 3 } y + A^(y) 

- y=V2/3 



14 



Note that tfi{Ai} s = A[ 1} (y) + (1/2) A^\y). The tunc- 

(2) 

tions An ' {(3) are given by 



71 7T /„ A? 



■dX , 



(A45) 



(2) 

and similarly for C„ (/?)■ We used the following iden- 
tity /* A n {a{?))? 2 df = (27r/3)Al 2) (/3), and similarly for 



Cn\f3). The explicit form of the functions A^ ! {0) and 
(/3) and their asymptotic formulas are given in • 
The asymptotic formulas for the tensors o§ and bfj k 
for a weak mean magnetic field B <C B cq /A are given by 

(SliVj + fyViXll (u 2 ) (°) + 3 <b 2 )(°)) 



45 



4/2 

ft n = — 5. 
0* 135 



8^(n-v)(3(u 2 )(°)-(b 2 )(°) 

(11 - e) fy<5 ifc - 2(2 - 7e) fi fc fy 



(A46) 



(A47) 



and for £> ^> B eq /4 they are given by 



ll/ 2 
3/3 

7^2 ' 3 
J 2 , 



5 y (n.A( fl ))(l-1.3c)-^ 



((u 2 ) + (b 2 )), 



ijk 



3 /3 



(1 - e) fij^fe + 5(1 - e) ClkSi 



&a (n-v) 



(A48) 



(A49) 



Using Eqs. f&25|l - l|A27|l and jA42 | -jA43 ]l we derive for- 
mulas for the contributions to the coefficients defining the 
nonlinear electromotive force due to a uniform rotation. 
In particular, the isotropic contribution to the hydrody- 
namic part of the a effect caused by a uniform rotation 
is given by 



o x 
a o. 



13 ) 



(A50) 



where a n is given by Eq. JTBJ, and the quenching func- 
tions 4^{B) and which determine af 1 , are given 
by 

$(B) = * 5 {Ax + A 2 - C x - C 3 }^ p , (A51) 
cf> 2 \B) = * 6 {Ax +A 2 -Cx- C 3 }^ p . (A52) 

The coefficients defining the nonlinear electromotive force 
due to a uniform rotation for a weak mean magnetic field 
B <C -B eq /4 are given by: 



f r 

45 



(aV,+^V l )(ll(u 2 )(°)+3(b 2 )(°)) 



;^(^-v)(3(u 2 )(°)-(b 2 )(°); 



(A53) 



21 



^(i-c)n 



Hit 

135 



1 13 
1 + Y e 



(A54) 
(A55) 



and for B ^> B cq /A they are given by 

«£■ « -^[in 2 /3(^A (s) )(l- £ ) 
+r 2 (^.V)((u 2 ) + (b 2 )) 



+^V l )((u 2 ) + (b 2 )), 

17 1a-)n, 



14^2/3 
i 2 



8-|(l-e)(r»^ lfe + r! fe ^) , 



(A56) 
(A57) 

(A58) 



r£ = 0(£1 2 ), and we took into account that (u 2 



(b 2 )(°) + Q(VB). Asymptotic formulas l|A4 



( u 2)(0) 

(|A47(I and (| A53|) - (| A55|) in the limit of a very small mean 
magnetic field coincide with those obtained in |2^| for 
q = 5/3. 

Now we study the effect of the mean differential rota- 
tion on the nonlinear electromotive force. We take into 
account the tensors If jmn (tJ), Ef jmn {\J) and Jg mn (U) 
in Eqs. HA1|) - (|A3|) . The contribution, Ef , to the nonlin- 
ear electromotive force caused by a mean velocity shear 
is determined by 



Jmnpqdpq + iT(k-B) [I mnpq fpq 



1 + 2^ 



r(h,a) 



. j(flV 



dk 



(A59) 



(for details, see |2(J), where the source terms — 



4(4), 4^ = 4^) and = i*(/; 



are 

determined by Eqs. (|A4|) - (|A6|) . in which fy, hij, are 
replaced by the corresponding correlation functions 
hfj , gfj that describe the contributions caused by a mean 
velocity shear. After the integration in Eq. I|A59|> . we 
obtain 



(A60) 



The tensor afj for an inhomogeneous turbulence is given 
by Eq. I|A67|I below. For a homogeneous turbulence afj = 
0. This case has been considered in j2||. The tensor b\ 
is given by 



\ 3 k 



7 

^ Qn $ijk 
n=X 



(A61) 



(see where the coefficient Q 3 = 0, and the other 

coefficients calculated for q = 5/3 are given by 



15 



Qi = 



Q2 



A? 3A? i8C« + e(>) + ^ + |cf>) + *i{.4i + 2A 2 + |d - ejW + A 2 + fc\ 



+* 2 ^ -Ax + -d + e{A t - 5Ci) }• - (1 - e)* 3 {Ci} - *o{2^i - 3d} 



(A« + 4 2) + 4Cf >) + e( + 4 2) + f Cf 



74 / 61 



64 - 



+* 2 {A 1 - 27Ci - e(Ai - 35Ci)} + (1 - e) \^ 3 {-2A 1 + 7d} - — * 4 {Ci} + 16* 5 {d} 
11 



-^ {2A 1 -—C 1 



Qa = 



Q 5 = 



(1-C) *2{Ai+Ci}-* 3 {Cl} 



4 2) +4 2) -y^ 2) + e(4 2) 



,(2) 



26 r (2) 



*o{Ci} 

*i<! Ai - A 2 



-Ci -2e Ai-i 



+-Ci ) [> + (!- e)( * 2 {Ai + d} - * 3 {d} ) + *o{Ci} 



,(2) 



AC 



(2) 



32 



,(2) 
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-3^ 2 + yd + 2eU 2 -yC 3 



(27- 35e)* 2 {C 3 } 



-f] , ; ( ^ 3 {A 2 - 7C 3 } + y* 4 {C 3 } - 16* 5 {C 3 }) + *o(a 2 - yC 3 



4 (2) 14 r (2) 

^2 y ° 3 



e 3A 



(2) 



yCf } ) +$i|A 2 + ^(14 < iC, (I -<)( >T' 2 {2.! 2 r :i } 



* 3 {A 2 + C 3 } + *o{^2 + C 3 } 



Here 

^ijfc = £ ijp{9U) p k , = Sikp(dU) p j 



? (3) 

?(6) 



e jkp (dU) 



j(7) 



The coefficients defining the shear-current effect and the 
nonhclical a effect are determined by 



Co 



2Q 4 + Qe + 2Q 7 ) 



1 



-co - +2Qs) 



(A62) 
(A63) 



Thus, the nonlinear coefficient cto(-B) and <Ji{B) are de- 
termined by 

a (B) = ^ a {A 1 + A 2 } + * 6 {Ci + C 3 } , (A64) 

fi(-B) = -<r (B) + *c{^i} + *«i{^2} + *e{Ci} , 

(A65) 

where 

*a{X} = \ 



(l + e)A( 2 )(V2/3) + [* -(l-e)(*i 



* C {A} 



*2 + * 3 )]{^} 

(3e - 13)X (2) (V2/3) + [12* 2 - 4* 

-16*i + (1 - e)(57*i - 51# 2 + 9* 3 
-32* 4 + 24* 5 )]{A} , 

-{l + e)xW(V2p) + y Q {X} 

2A^(V2/?)-3(l-e)*i{X} 

(34 + 12e)X (2) (V2/3) + [4* 2 - 6* 

20#i - (1 - e)(*i + 9*2 - 3* 3 )]{A} 



and the functions *fc{A} are given by 

*o{A} = - l -{l + e)X^(Q) + {2-e)X^{V2f3) 



16 



-_(l- e )X(2/3 2 ), 



Z7T 



V 2 {X} = ZX^\V2p) 
* 3 {X} = -6X<- 2 \V2/3) 



2tt 



X(y) + -yX'(y) 



y=2f3 2 



2tt 



yX\y) 



y=20 2 



AX (2 \V2(3) - - 

1 



2X(y)+yX'(y) 



y=2(3 2 



* 5 {X} = 



X(y) 



+ -yX l (y)+y 2 X"(y) 
The tensor is given by 



(A66) 



/2 



Fi «5<j (W-AW) + F 2 WiA^ + F 3 W 3 A 



+F 4 Sfl A£> + F 5 S^A^ + e[F 6 5 i3 - (W-A 



+F 7 W t kf + F a WjA^ + F 9 s\f n A^ 
+F m sgA^} , 

where 

Fx = -HW){A 1 } + ~[hW{A 2 } 

+A (2 \V2P)}, 
F 2 = -i[(6GW-3^){A!}-4 2) (V2/3)] 

F 3 = -l[ff<2>{Ai+A 2 }+A< 2 >(V20) 
+A¥\V2(3)], F 4 = -2F 



(A67) 



-2F, 



3 j Fg — — 2Fs 
(3G (2) +8ff (2) ){Gi + G 3 } 



+4[cf ) (V2/?)+^ 2) (^/3)] , 

F 6 = -3(G( 2 )-ff( 2 )){^}-i[^( 2 ){A 2 } 

-4 2) (y2/3)], 
F 7 = i[(6G( 2 )-7 J ff( 2 )){A 1 } + 4 2) (V2/3)] , 

F 8 = ±IH^{A 1+ A 2 }~a{ 2 \V2P)-A¥\V2P)} 
F w = -2F 6 - 4(G (2) + 2H^){d + C 3 } 
+ y[Cp ) (v/2/3)+cf (V2/3)] , 



and 

G (2) {X} - 10X( 2 )(V2/3)-^-[6X(y) 

47T 

H^{X} = 4A^(\/2/3) - — X{2(3 2 ) 

2tt 

= ^{X}^ , 

G m {X ) - H W{x} = * 3 {x}^ , 

For the derivation of Eq. 1A67(I we used the following 
identities 

£inm,A n I3jKj m pqX7pUq — 2(Gi ~t~ Cz)^-inqA n {^dU^jqjBj , 
{SimqKjmSpn — Ei n qKjp)A n BjSlpUq — (Ai + A 2 )ApBj 

x e ijq (dU) pq - e ipq (dU)j q - \{k P W + 8 i3 W p ) , 



SijWn) 



A n Bj . 



An additional contribution to the isotropic part (pn 3 oc 
aSij ) of the nonlinear a effect [see Eq. Ijl5(l ] due to both, 
inhomogeneity of turbulence and mean differential rota- 
tion in a nondimcnsional form in spherical coordinates is 
given by 



where 



# 6 {Jf }A (w) + e^ 7 {X}A^ 



xsm0—(6Q) 



X=Ci+C 3 

(A68) 



* 6 {A} = 



* 7 {X} 



-f^ (2) ^) + li 



34A(y) 



-yx'(y) 



10X(y) 



-yx'(y) 



The contribution to the nonlinear a effect due to both, 
inhomogeneity of turbulence and mean differential rota- 
tion for a weak mean magnetic field B <C F eq /4 is given 
by 



a&n = 2LW* 



9 L„ 



3 



sin O-^m , 



(A69) 



and for B ^S> B eq /4 it is given by 



JQ 



1 LW» 



= -7^9 -T^ eA( ' } sine— . (A70) 



9/3 



86 
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Equations for a in cylindrical coordinates can be Note that the a term has been also calculated in |4fJ 
obtained from Eqs. i|A68[) - HA70|) after the change for a kinematic problem using the second-order correla- 
sa\9{d/d9) — > p{d/dp). tion approximation (SOCA). 
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